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^ : ABSTRACT 
(N 

CN ■ This paper is aimed at giving an overview of the global properties of the 



rich cluster of galaxies ABCG 209. This is achieved by complementing the 

, already available data with new medium resolution spectroscopy and NIR 

OO . 

■ photometry which allow us to i) analyse in detail the cluster dynamics, dis- 



tinguishing among galaxies belonging to different substructures and deriving 
' their individual velocity distributions, using a total sample of 148 galaxies in 

5t 1 the cluster region, of which 134 belonging to the cluster; ii) derive the cluster 

NIR luminosity function; iii) study the Kormendy relation and the photomet- 
ric plane of cluster early- type galaxies (ETGs). Finally we provide an extensive 
photometric (optical and NIR) and spectroscopic dataset for such a complex 
system to be used in further analyses investigating the nature, formation and 
evolution of rich clusters of galaxies. 

The present study shows that the cluster is characterised by a very high 
value of the line-of-sight velocity dispersion: CTv = 1268^81 km s^i, that results 
in a virial mass of M^ir = 2.95^0 78 x 10^^ Kj^ within R^ir = 3.42/i^g-^ 
Mpc. The analysis of the velocity dispersion profile shows that such high 
value of (T„ is already reached in the central cluster region. There is evidence 



* Based on observations collected at European Southern Observatory (ESO ObsID.s 068.A-0116, 074.A-0073 and 075.A-0845) 
and at the Telescopio Nazionale GaUleo (TNG ObsID AOT12/TAC.07). 
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of three significant substructures, the primary one having a velocity dispersion 
of CTv = 847j^4g km s"^, which makes it consistent with mass estimates from 
weak lensing analyses. 

This observational scenario confirms that ABCG 209 is presently under- 
going strong dynamical evolution with the merging of two or more subclumps. 
This interpre tation is also supported by the detection of a radio halo (Gio- 



vannini et al. 



20061) suggesting that there is a recent or ongoing merging. 



Cluster ETGs follow a Kormendy relation whose slope is consistent with 
previous studies both at optical and NIR wavelengths. We investigate the 
origin of the intrinsic scatter of the photometric plane due to trends of stel- 
lar populations, using line indices as indicators of age, metallicity and a/Fe 
enhancement. We find that the chemical evolution of galaxies could be re- 
sponsible for the intrinsic dispersion of the Photometric Plane. 

Key words: galaxies: clusters: general - galaxies: clusters: individual: ABCG 209 
- galaxies: distances and redshifts - galaxies: kinematics and dynamics - galax- 
ies: photometry - galaxies: luminosity function 



1 INTRODUCTION 



Galaxy clusters are complex systems involving a variety of interacting components: galaxies, 
ho t and cold gas, dark matter. Among them, about 50% are unrelaxe d syst ems (Smith et 
2OO5I ) and a large fraction contain substructures (e.g., Girardi et al. 



al. 



2006 and references 



therein), suggesting that they are evolving via merging processes from poor groups to rich 
structures. In fact, in hierarchical clustering cosmological scenarios, galaxy clusters form 
from the accretion of subunits. Numerical simulations show that clusters form preferentially 
through the a nisotr opic accretion of subclusters along filaments (e.g., Colberg et al. Il999l . 



Diaferio et al. 



200l|). 



Detailed multi-band studies of such systems are crucial to probe structure formation 
scenarios as well as to investigate galaxy evolution processes. ABCG 20 9 is a n ideal target 
for such a study, since it is a rich (richness class R = 3; A bell et al. Il989l ). X-ray lumi- 
no us (Ly fO. 1-2.4 keV) ~ 2.7 x 10^^ /ifo^ ergs^S Ebeling et al. 



al. 



19961 : Tx ~ 10 keV, Rizzaet 



19981 ) ■ and massive c 



Paulin-Henriksson et al. 



uster (M(R < R^ir) = 2.3-3.1 x lO^^/ifo^M© Mercurio et al. 



2003al . 



20071 ). Evidence for a complex dynamical status comes from the 



t Tables 2 and 3 are available only in electronic format. 
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X-ray emission: the hot g as is elongated and distributed asymmetrically showing two main 
clumps (Rizza et al. Il998l ). However no strong cooling flow is detected. Moreover, the young 
dyna mical state is indicated by the possible presence of a radio halo (Giovannini et al. 



19991 . 



20061 ). which has been suggested to be the result of a recent cl uster r nerge r, through the 



acceleration of relativistic particles by the merger shocks (Feretti 



2002 



20071). 



uster w ere analysed through a spectroscopic survey of 112 



Dynamical properties of the c 
cluster members (Mercurio et al. 
were derived out to radii of 3- 4 h^n Mpc, taking advantage of wide field C 



2003aD and photometric properties of the cluster galaxies 



i?-band images (E 
La Barbera et al. 



aines et al. 



2003a 



20041). Our analyses (Merc urio et al. 



2003 



2003a 



^HT B- and 



2003b . 



2004 



2Qm showed that ABCG 209 has a 



20041 : Haines et al. 

very complex structure, with: a) a high value of the line-of-sight velocity dispersion, with 
(Jt, = 1394^99 km s^^; b) a significantly non-Gaussian redshift distribution, and a velocity 
gradient (and elongation of the cD galaxy) in the SE-NW direction; c) the presence of 
substructures in the X-ray emission; d) the presence of a merger clump observed 1-2 Gyr 
after the merging; e) a strong spatial and spectral segregation of galaxies, with (i) young 
blue emission line galaxies uniformly distributed in the low density regions, (ii) blue post- 
starburst galaxies aligned in a direction perpendicular to the cluster elongation, and (iii) 
red post-starburst galaxies and early type galaxies distributed along the cluster elongation. 
Moreover, we have examined the effect of cluster environment, as measured in terms of the 
local surface density of R <23.0 mag galaxies, on the global properties of the cluster galaxies, 
through their luminosity functions, colour-magnitude relations, and average colours. The LFs 
for galaxies within the virialized region are found to be well described by single Schechter 
functions, although there is an indication of a dip at i?=20-20.5 mag. The faint-end slope 
shows a strong dependence on environment, becoming steeper at more than 3cr significance 
level from high- to low-density environments. We explain this trend as a combination of the 
morphology-density relation and dwarf galaxies being cannibalised and disrupted by the cD 
galaxy and interactions with the intra-cluster medium i n the cluster core. 

Our weak lensing analysis (Paulin-Henriksson et al. 120071 ) confirms that ABCG 209 is a 
massive cluster, although the mass estimated by weak len sing (M (R < R200) = 7.7l;|;7 x 
IO^^/i^o^Mq) is lower than that obtained by Mercurio et al. ( l2003al ) from the analysis of the 
dynamical properties. The centres of the X-ray emission, dark matter and galaxy distribu- 
tions all appear offset from one another, with the centre of mass found from the weak lensing 
analysis lying between that of the X-ray and galaxy distributions, with all three centres of 
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mass aligned on the main SE-NW axis of the cluster (see Figs. 6 and 7 of Paulin-Henriksson 



et al. 



20071 ). Such an effect is seen for the more extreme "Bullet cluster" (Clowe et al. l2004l ). 



and seems to reflect the different responses of the gas and dark matter components to the 
merger, in agreement with the merging scenario for ABCG 209. 

The above results show that it is crucial to relate the properties of member galaxies to the 
global properties of clusters, such as mass and dynamical state. The spatial and kinematical 
analysis of member galaxies allows to detect and measure the amount of substructure and 
to identify and analyse possible pre-merging clumps or merger remnants. In addition, sub- 
clustering is important for processes of galaxy evolution. In fact, although the high galaxy 
velocity dispersion associated with a relaxed, virialized clusters inhibits galaxy mergers, in 
subclusters, where the velocity dispersions are lower, such mergers a re mo re probable and 
could have an important role in building elliptical galaxies (see Moss l2006l ). 

Identifying substructures is difficult with galaxy positions alone, since projected distri- 
bution of galaxy clusters could show clumps that are not real. The optical spectroscopy of 
member galaxies is the most reliable tool to investigate these clumps and allows us to study 
the dynamics of cluster mergers, since it provides direct information on the velocity field. 
However, this is often an arduous investigation due to the limited number of galaxies usually 
available to trace the velocity field. 

In order to achieve the resolution needed to understand the complex dynamics of ABCG 209 
and to investigate the environmental effects on galaxy evolution, we complement the pre- 
viously described dataset with medium resolution spectra analysing a total sample of 148 
galaxies in the cluster region. The enlarged data sets allow us to better distinguish among 
galaxies belonging to different substructures and to derive their individual velocity distribu- 
tions. A detailed study of substructures is also important to investigate the reason for the 
discrepancy between lensing and dynamical mass estimates. We also present new K-band 
observations, that allow us to derive the cluster NIR luminosity function (LF) and are di- 
rectly related to the galaxy mass function. Additionally, at infrared wavelengths the galaxy 
luminosities do n ot de pend strongly on the details of their stellar populations (Gavazzi, 



Pierini, & Boselli 



19961). 



In this paper we provide the spectroscopic catalogue of the new targets and a photo- 
metric catalogue with measured photometric redshifts, V-, R- and i^"-band total mag- 
nitudes and spectroscopic redshifts when available. We also derive structural parameters of 
cluster galaxies, selected according to the photometric and spectroscopic redshifts; and we 
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investigate the photometric correlations of ETGs, such as the rnean surface brightness-size 
relation, also known as the kormendy relation (KR; Kormendy 119771 ) and the Photometric 
Plane (PHP; La Barbera et al. boosi). 

The paper is organised as follows. We present the new data in Sect. 2 and we derive 
photometric and spectroscopic redshifts in Sect. 3 where the catalogue of the cluster members 
is also described. We analyse the global cluster properties in Sect. 4 and the internal dynamics 
and substructures in Sect. 5. We describe the derivation of structural properties of ETGs in 
Sect. 6, deriving also the KR and the PHP. Finally we discuss and summarise our results in 
Sect. 7. 

Unless otherwise stated, we give errors at the 68% confidence level (hereafter c.l.). 
Throughout the paper, we assume a flat cosmology with fim = 0.3, = 0.7, and Ho=70 km 
s~^ Mpc~^. For this cosmological model, 1' corresponds to 271 kpc at the cluster redshift. 



2 OBSERVATIONS AND DATA REDUCTION 

Spectroscopic observations were carried out at the ESQ New Technology Telescope (NTT) 
with the ESQ Multi Mode Instrument (EMMI) and at the Telescopio Nazionale Galileo 
(TNG) with the Device Optimized for the LOw RESolution (DOLORES), while NIR pho- 
tometric data were collected with the Son OF ISAAC (SOFI) at NTT. 



2.1 Photometry 

New fT-band images for the cluster of galaxies ABCG 209 were collected with SOFI oper- 
ating in the LARGE FIELD observing mode, providing a field of view of 5' x 5' with a pixel 
scale of 0.288"/pxl. Eight overlapping fields were observed, each of them with a dithered 
sequence of 36 exposures. We adopted a 20" dithering box, with DIT and NDIT values of 
6 s and 10, respectively. This resulted in an integration time of 60s for each exposure, and 
a total integration time of 2160s for each field. Sky conditions were photometric during five 
o ut of six observing n ights. During those nights we observed standard stars from the list 



of 



Persson et al 



(119981 ) ■ at five different positions on the chip. 
The data were reduced using FORTRAN routines developed by the authors. For each 
dithering sequence, the exposures were dark subtracted and flat-field corrected, using a 
superfiat frame obtained by median combining all the images taken during a given night. 
After this procedure, the magnitudes of the standard stars showed a rms variation across 
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the chip of ~ 0.05 mag. To achieve a better accuracy, we retrieved K-band images from 
the ESO archive for a standard star observed at different positions in a 5 x 5 grid on the 
SOFT frame. We chose standard star observations as close as possible in time to those of 
ABCG209. Illumination correction frames were then obtained for each night by measuring 
the magnitude of the standard star as a function of the position on the frame. This allowed 
the low spatial frequency component of the flat-field to be corrected to better than 1%. Since 
sky subtraction is a very troublesome step for the reduction of NIR data, particularly in 
high density regions, this point was carefully dealt with by a two step procedure. First, each 
exposure was sky subtracted by computing the sky frame from the median of the six closest 
frames along the sequence. The images were then registered with integer shifts and combined 
using a sigma clipping algorithm for cosmic ray rejection. This procedure alone overestimates 
the sky level in the extended halos of galaxies. To minimise this effect, we reiterated the 
sky subtraction as follows. The initial combined images were used to obtain mask frames 
for the sources in the field, by running SExtractor with the checkimage OBJECTS option. 
For each sequence, the mask was expanded in order to 'cover' the galaxy halos and was 
de-registered to each dithered exposure. The hot and bad pixels were also masked. The sky 
frames were estimated as the average of the six closest frames to the exposures from that 
sequence, rejecting masked pixels. The exposures of each dithering sequence were then sky 
subtracted and combined with the IRAfQ task IMCOMBINE. The images of the different fields 
were combined by taking into account their different zero-points, resulting in a final i^-band 
mosaic with an average seeing FWHM of ~ 0.8". The mosaic is shown in Fig.^ 



The photometric calibration was performed for the Kg standard filter (jPersson et al 



19981 ). deriving the instrumental magnitudes of the standard stars within an aperture of 



8" diameter. The airmass correction was performed by using an extinction coefficient Ak 
= 0.05 mag/airmass, which was derived by comparing the magnitudes of bright objects in 
the field for different airmasses. The typical accuracy on the zero-point of each photometric 
night amounts to ~ 0.01 mag. The images of the different fields were normalised to the 
average zero-point value of 22.425 (scaled to 1 s exposure time). 



^ IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Figure 1. K-band mosaic image of ABCG 209. Contours arc obtained from tlic projected number density map of K < 17.5 
mag galaxies. The contour density values are equally spaced by 1.27 galaxies arcmin"^, with the outermost and innermost 
contours corresponding to surface density values of 1.4 galaxies arcmin"^ and 9 galaxies arcmin"^, respectively. The mosaic 
covers an area of ~ 0.0525 square degrees, corresponding to ~ 1.39 h^g'^ Mpc'^ at z=0.209. 



2.2 Spectroscopy 

The spectroscopic data were obtained with the multi-object spectroscopy (MOS) mode. 
Targets were selected according to the B-R vs. R colour- magnitude relation (CM). We 
gave priority to the galaxies lying on the CM relation up to i?=20.0 mag. We acquired 
one mask with EMMI (field of view 5' x 8.6'), centred on aj20oo =01:31:53.0, 5j2ooo =- 
13:37:42.8, allocating 28 slits and one with DOLORES (field of view 6' x 7.7'), centred 
on q;j2ooo =01:31:52.0, (5j2ooo =-13:37:30.8, allocating 33 slits. We integrated for 8100 sec 
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with EMMI-Grism#5 yielding a dispersion of ~1.6 A/pxl (resolution=~4.8 A FWHM) in 
the spectral range 380-702 nm, and for 16200 sec with DOLORES-grism MR-B, giving a 
dispersion ~1.9 A/pxl (resolution=~4. 2 A FWHM) in the spectral range 350-700 nm. 

Each scientific exposure (as well as flat fields and calibration lamps) was bias subtracted. 
The individual spectra were extracted and flat field corrected. Cosmic rays were rejected in 
two steps. First, we removed the cosmic rays lying close to the objects by interpolation 
between adjacent pixels, then we combined the different exposures by using the IRAF task 
IMCOMBINE with the algorithm CRREJECT (the positions of the objects in different exposures 
were checked before). Wavelength calibration was obtained using He-Ar and He lamp spec- 
tra for EMMI and DOLORES, respectively. The typical rms scatter around the dispersion 
relation was ~ 8 km s^^ and 11 km s~^ for EMMI and DOLORES respectively. 

The positions of the objects in the slits were defined interactively using the iRAF package 
APEXTRACT. The exact object position within the slit was traced in the dispersion direction 
and fitted with a low order polynomial to allow for atmospheric refraction. The spectra 
were then sky subtracted and the rows containing the object were averaged to produce the 
one-dimensional spectra. The signal-to-noise ratio per pixel of the one-dimensional galaxy 
spectra belonging to the cluster ranges from about 6 to 18 in the region around the Mg 
line (at 625 nm in the spectra at z~0.2) for EMMI spectra and from about 8 to 27 for the 
spectra obtained with DOLORES. 



3 REDSHIFT MEASUREMENTS 
3.1 Spectroscopic redshifts 



Redshifts were derived using the cross-correlation technique (ITonry &: David Il98lf). as im- 



pleme nted in the RVSAO package. We adopted galaxy spectral templates from 



Kennicutt 



(119921 ). corresponding to morphological types E, SO, Sa, Sb, Sc and Ir. The correlation was 
computed in the Fourier domain. We define the redshift as the value given by the best tem- 
plate spectra, i.e. the template producing the highest value of the correlation R — value given 
by RVSAO that gives an indication of the signal-to-noise of the correlation peak. 

Of the 61 observed spectra, 44 turned out to be at the redshift of ABCG209, 1 is a 
foreground galaxy and 6 are background galaxies. In 10 cases we could not determine the 
redshift. In Table [H we report spectroscopic redshift measurements of the observed galaxies. 

In order to estimate the uncertainties in the redshift measurements, we considered the 
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Table 1. Spectroscopic data. Running number for galaxies in the present sample, right ascension and declination (Col. 1, 
Col. 2); heliocentric corrected redshift z (Col. 3). 



01 31 
01 31 
01 31 
01 31 
01 31 
01 32 
01 31 
01 31 
01 32 
01 31 
01 31 
01 31 
01 31 
01 31 
01 31 
01 31 
01 31 
01 31 
01 31 
01 31 
01 32 
01 31 
01 31 
01 31 
01 31 
01 31 
01 31 
01 31 
01 31 



45.33 
56.81 
47.02 
48.41 
45.78 
01.00 
52.41 
42.47 
01.27 
57.49 
47.06 
45.35 
46.73 
52.37 
52.91 
56.40 
51.79 
53.25 
52.52 
42.34 
04.29 
49.42 
44.94 
49.39 
49.14 
48.95 
48.33 
43.11 
40.63 



-13 41 53.0 
-13 41 31.7 
-13 41 17.9 
-13 41 00.4 
-13 40 37.8 
-13 40 17.5 
-13 40 02.0 
-13 39 44.3 
-13 39 33.3 
-13 39 23.4 
-13 39 06.7 
-13 38 46.0 
-13 38 29.2 
-13 37 40.6 
-13 37 34.6 
-13 37 18.3 
-13 36 54.8 
-13 36 44.2 
-13 36 27.5 
-13 36 06.4 
-13 35 57.0 
-13 35 44.5 
-13 35 11.3 
-13 34 59.3 
-13 34 45.5 
-13 34 34.9 
-13 34 23.5 
-13 34 07.5 
-13 33 46.3 



0.2004 ± 0.0004 
0.2182 ±0.0002 
0.2030 ± 0.0002 
0.2079 ± 0.0002 
0.2031 ± 0.0002 
0.2113 ±0.0003 
0.1989 ±0.0003 

0.2539 ± 0.0003" 
0.2103 ±0.0005 
0.2165 ±0.0003 

0.3116 ± 0.0003" 

0.2522 ± 0.0003" 

0.2519 ± 0.0003" 
0.1980 ±0.0002 
0.2150 ±0.0002 
0.2080 ± 0.0003 
0.2074 ± 0.0002 
0.2098 ± 0.0004 
0.2110 ±0.0003 
0.2080 ± 0.0003 

0.3702 ± 0.0005" 
0.2164 ± 0.0004 
0.2130 ±0.0002 
0.2078 ± 0.0003 

0.4025 ± 0.0002" 
0.2058 ± 0.0005 

0.1486 ±0.0004^ 
0.2109 ±0.0002 
0.2097 ± 0.0003 



^ Background galaxy. 
^ Foreground galaxy. 

error calculated with the cross-correlation technique, which is based on the width of the 
peak and on the amplitude of the antisymmetric noise from the cross correlation (cf. manual 
of the XCSAO task). The wavelength calibration errors (see Sect 12. ip turned out to be negli- 
gible in this respect. The errors derived f rom t he cross-correlation could however be smaller 
than the true errors (e.g. Boschin et al. |2006| and references therein). We checked the er- 
ror estimates by comparing redshifts computed for 17 galaxies observed with EMMI and 
5 with DOLORES which were already acquired in a previous observation carried out with 



EMMI ( iMercurio et al.ll2003al ). The two data sets agree with a one-to-one relation, and a 
reasonable value of for the fit, in particular for the spectra acquired with DOLORES, was 
obtained when the errors derived from the cross-correlation were inultiplied by a correction 
factor ~1.5. A similar correction w as obt ained by Malumuth et al. (119921 : 1.6), Bardelli et al. 
(1994; 1.87), and Quintana et al. f 200cj : 1.57). Since we are comparing redshifts measured 
from spectra acquired in different periods and with different instruments, the multiplicative 
factor obtained through this comparison takes also into account possible external errors. 
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3.2 Photometric redshifts 



The redshifts of galaxies without spectroscopic information were estimated 
photo metri c redshift technique, usin g a similar approach to that described in 



3V means of the 



Busarello et al. 



fl2002h and 



La Barbera et al 



(l2003d ) . Photometric redshifts were estimated according to the 



Massarotti et al. 



2001 



and references 



Spectral Energy Distribution fitting method (see 
therein), using the new i^'-band photometry, as well as the already published B-, and 



2004, 



i?-band photometry of ABCG 209. The B- and -R-band data, described in Haines et al 
are available on the entire f ield cov ered by the new i^-band photometry, while the V^-band 
imaging (see Mercurio et al. l2003bl ) covers a smaller region of ~ 160 arcmin^ centred around 
the cluster core. Photometry of so urces in the i^"-band mosaic are obtained by means of 
SExtractor (Bertin & Arnouts ll996l ). as detailed in Sec. 14.1.11 Galaxy colours were measured 
within an aperture of 5" diameter. Since all the BVRK images have similar seeing FWHM 
values (0.7" — 0.9"), no aperture corrections were applied to galaxy colours. 

Photometric redshifts were measured for the 642 galaxies, with available photometry 
in B, R, and K bands, brighter than R = 21 mag: the magnitude limit for spectroscopic 
observation. We looked for redshifts in the range z G [0.0, 1.0] with a step of 0.01, impos- 
ing that at a given redshift galaxy templates would be younger than the age of the uni 



verse in the adopted cosmology. We used the GISSEL03 spectral code (IBruzual fc Chariot 



20031 ) to produce galaxy templates with a Scalo IMF and an exponential SFR, e The 



colours of E/SO, Sa/Sb, and Sc/Sd spectra were modelled by choosing r = 1,4 and 15 Gyr, 
respectively, while early-type galaxies with different metallicities were described by us- 
ing E/SO models with Z/Zq=0.2, 0.4, 1 and 2.5. The differential dust extinction of the 
Milky Way w as included in the computati o n of model colours by adopting the extinc- 



tion curve of 



Cardelli. Clavton and Mathia (1l989l ) and a colour excess of E(B — V) 



0.019 (jSchlegel. Finkbeiner and David Il998h . The uncertainty on the photometric redshift, 
6z, was estimated by performing numerical simulations, shifting galaxy colours according to 
their corresponding uncertainties, and recomputing each time the photometric redshifts. 

We compared photometric and spectroscopic redshifts for = 134 galaxies in the cluster 
field. We found that galaxies with large relative uncertainties on photometric redshifts were 
biased toward higher redshift values. In order to reduce this systematic effect, we removed 
from the sample those galaxies with a value of 6z/l + z larger than 0.07. This selection leads 
to a sample of 399 (out of 642) galaxies with reliable photometric redshift estimates (see 
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Figure 2. Distribution of photometric redshifts for galaxies with a relative uncertainty on photometric redshifts smaller than 
0.07. The dashed lines mark the redshift range that defines likely cluster members. 



Table [2]). As shown in Fig. [2|, the distribution of these photometric redshifts is dominated 
by the peak around z ~ 0.2, indicating that most galaxies with i? ^ 21 in the i^-band field 
are actually cluster members. After comparison with spectroscopic redshifts, we chose to 
mark as likely cluster members those galaxies in the photometric redshift range of 0.14 to 
0.26. This selection leads to a sample of 292 (out of 399) galaxies which are likely cluster 
members. 



3.3 Catalogues 

The new spectroscopic catalogue is presented in Table [H which contains the sky coordinates 
a (Col. 1) and 5 (Col. 2), and the spectroscopic redshift (Col. 3). 

The catalogue of galaxies with reliable photometric redshifts { 5z/l + z < 0.07) is pre- 
sented in Table 121 available in electronic format, which includes the identification number of 
each galaxy, ID (Col. 1), right ascension and declination (J2000), a and 5 (Cols. 2 and 3), 
B (Col. 4), V (Col. 5) and R magnitudes when available (Col. 6), K magnitudes (Col. 7), 
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Table 2. Photometric data. Running number for galaxies in the present sample, ID (Col. 1); right ascension and declination 
(Col. 2, Col. 3) B-, V-,R- a^id K-band Kron magnitude (Col. 4, Col. 5, Col. 6, Col. 7); heliocentric corrected redshift z (Col. 8); 
photometric redshift (Col. 7). 



ID 




a 




s 






B 




V 


R 


K 


Zgpec 


^phot 


1 


01 


31 


47.23 


-13 


35 


31.2 


21 


97 ± 0.02 


20 


35 ±0.03 


19.73 ±0.01 


17.06 ±0.02 




0.155 


2 


01 


31 


49.46 


-13 


37 


26.9 


19 


12 ± 0.01 


18 


02 ±0.02 


17.57 ±0.01 


15.11 ±0.01 


0.2140 ± 0.0002 


0.170 


3 


01 


31 


55.10 


-13 


37 


24.1 


21 


92 ± 0.01 


21 


27 ±0.03 


20.63 ±0.01 


18.71 ±0.09 




0.451 


4 


01 


31 


55.68 


-13 


37 


19.9 


23 


09 ± 0.03 


21 


43 ± 0.03 


20.94 ±0.01 


18.18 ±0.07 




0.147 


5 


01 


31 


51.56 


-13 


37 


18.0 


22 


03 ± 0.02 


20 


46 ± 0.02 


19.67 ±0.01 


16.54 ± 0.01 




0.209 


6 


01 


31 


51.28 


-13 


37 


17.5 


22 


54 ± 0.02 


21 


06 ± 0.02 


20.34 ± 0.01 


17.47 ± 0.03 




0.160 


7 


01 


31 


56.40 


-13 


37 


18.2 


22 


17 ±0.02 


20 


65 ± 0.02 


19.97 ±0.01 


16.93 ± 0.02 
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heliocentric corrected spectroscopic redshift, z^pec (Col. 8), if available, and the photometric 
redshift, Zphot (Col. 9). 



4 GLOBAL CLUSTER PROPERTIES 

4.1 NIR luminosity function of cluster galaxies 

4.1.1 K-hand Magnitudes and Completeness 



To perform source detection, we ran SExtractor (Bertin & Arnouts Il996l ) on the i^'-band 
mosaic of ABCG 209, obtaining a i^-band catalogue with a total of 2628 objects. For each 
source, we measured aperture magnitudes within a 5" diameter aperture, and Kron mag- 
nitudes within an aperture of diameter a • tk, where tk is the Kron radius (Kron 1980). 
We chose a = 2.2, for which the Kron magnitude is expected to enclose 92% of the total 
flux, and we computed the total magnitude, K, by subtracting 0.08 mag from the Kron 
magnitudes. 

As shown in Fig. [3] (left panel), objects were classified as stars and galaxies according 
to the distribution of sources in the SI versus Kron magnitude diagram, where SI is the 
stellarity index parameter of SExtractor. We classified as stars those objects whose SI value 
is larger than a given threshold, S/min- The value of SImin was chosen by adding simulated 
stars and galaxies to the K-band mosaic, and measuring their SI and mx parameters by 
means of SExtractor, in the same way as for real sources. Simulated stars and galaxies were 
randomly generated in a magnitude range of = 14 mag to = 21 mag. In order to 
create the simulated images, we modelled real sources in the K-hand mosaic as the sum 
of three two-dimensional Moffat functions. The Moffat models were scaled according to the 
magnitude of simulated objects and added randomly to the i^-band mosaic. In the case of 
stars, the modelling was performed by fitting non-saturated sources with K < 16 mag and 
SI > 0.9. For galaxies, we fitted objects with SI < 0.2, in a magnitude range of -ft^min to 
Kma.x- We chose -ft'min = 17 mag and -ft'max = 18.5 mag, and we also varied the values of -ft'min 
and -ft'max in order to test the robustness of our completeness estimates (see below) with 
respect to these parameters. Fig. [3] shows the distribution of simulated objects in the SI 
versus K diagram. At magnitudes brighter than K ~ 17.5 mag, stars and galaxies are clearly 
well separated, with this separation disappearing at fainter magnitudes. In fact, ioi K ^ 17.5 
mag, simulated stars have a stellarity index ^ 0.4, while almost all simulated galaxies have 
SI <^ 0.8. In order to minimise the fraction of misclassified stars and galaxies, we adopted 
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Figure 3. Left panel: Stellarity index versus /sT-band magnitude diagram for all sources detected in the /sT-band mosaic of 
ABCG209 (red circles), for simulated galaxies (blue circles), and simulated stars (black dots). Shaded area indicates the 
selection criteria used to identify stars. Right panel: /f-band completeness function. Crosses, squares and triangles denote the 
completeness percentage as estimated by simulating galaxies in different magnitude bins (see the text) as shown in the lower 
left corner of the panel. The completeness percentage of stars is plotted by the circles. The solid curves represent the fits of the 
completeness functions by Fermi-Dirac functions. The dashed lines mark a completeness level of 90% for the galaxy catalogue. 

a value of S'/min that depends on the source magnitude, with S'/min = 0.95 for K ^ 17.5 
mag and Sl^i^ = 0.8 for K > 17.5 mag. The corresponding star/galaxy classification leads 
to a sample of 2334 galaxies in the i^-band field of ABCG209. The completeness of the 
i^-band catalogue was estimated by measuring the percentage of simulated galaxies and 
stars which are recovered by SExtractor as a function of K. The completeness functions are 
shown in Fig. [3] (right panel) for both stars and galaxies. The figure shows that the i^-band 
galaxy catalogue of ABCG 209 is complete at more than 90% down to K ~ 19 mag, with 
the completeness level falling rapidly to ~ 50% at 19.5, and to ~ 20% at ~ 20 mag. 
As one would expect, the stellar catalogue is characterised by a higher completeness at a 
given magnitude, being more than 90% complete down to ~ 19.8 mag. The figure also 
shows that the galaxy completeness function does not depend significantly on the magnitude 
range of the sources used to create simulated galaxies. In fact, changing the values of -ft'min 
and -ft'max, as shown in the figure, we found that the completeness estimates do not change 



significantly. 



4-1-2 K-band Luminosity function 

The K-band LF of ABCG 209 was derived in a square region of 4.7 arcmin side size around 
the cluster center, by computing number counts of galaxies brighter than K = 19.5 mag. 
Using the simulation results discussed in Sec. I4.1.H we corrected number counts for (i) 
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incompleteness, dividing by the completeness function shown in Fig. [3] (right panel), and 
(ii) contamination from misclassified stars and galaxies, dividing number counts by the ex- 
pected fraction of misclassified galaxies and subtracting the expected number of misclassified 
stars in each magnitude bin. Number counts were then corrected for background/foreground 
gal axies, by usiiig field galaxy counts from the Calar Alto Deep Imaging Survey (CADIS, 



see 



Huang et al.ll200ll ). a medium deep i^-band survey with a total area of 0.2 deg^ and 



a completeness magnitude of 19.75 mag. Since our i^-band photometry is 50% complete 
down to mx = 19.5 mag, the CADIS data are a suitable dataset for the estimation of field 
galaxy counts. The resulting LF of ABCG 209 and its fit with a single Schechter func- 
tion, are shown in Fig. |H The error bars on number counts take into account Poissonian 
uncertainties on both field and cluster counts. The LF fit was performed by a minimi- 
sation routine, accounting for the finite size of magnitude bins by integrating the Schechter 
function in each magnitude bin. Since brightest cluster galaxies usu ally prevent a good fit 



of t he LF to be achieved with a single Schechter function (see e.g. Ide Propris et al. 



and 



1998 



de Propris et al.lll999l ). the cD galaxy of ABCG 209 was excluded from the fitting. The 



Schechter function fit, whose parameters are the faint end slope, a, and the characteristic 
magnitude, K* , provides a good description of cluster galaxy counts. The only magnitude 
bin showing a larger deviation is that aX K = 17.5 mag, where a marginal (~ 1.5cr signifi- 



cant) deficiency of cluster galaxies is found. Using an R — K colour o: 
of early-type cluster galaxies at z ~ 0.2 (e.g. 



La Barbera et al. 



3.1, which is typical 



mm . hereafter LM04), the 



magnitude of = 17.5 mag can be transformed to i? ~ 2 0.6 ma g, that corresponds to the 
magnitude bin (20 < R <21 mag) where Mercurio et al. (l2003bl ) find some indication of a 
dip in the optical luminosity functions of ABCG 209. The best-fitting value of the ii'-band 
Schechter function parameters turn out to be a = —0.98 ± 0.15 and K* = 14.80 ± 0.30 mag, 
respectively. The uncertainties on a and K* were derived by randomly shifting galaxy num- 
ber counts according to their uncertainties, and then re-computing the best-fitting Schechter 
function. Joint probability contours for a and K* are shown in Fig. [3l The above value of 
a is fully consistent with those found from previous studies of NIR lu minosity functions 
of clu ster galaxies, such as a =~ 1.0 fide Propris et al. 
20001), and a ~ -1.2 (LM04). 



19981), a - -1.3 flAndreon fc Pello 



16 A. Mercurio et al. 



-1 1 1 r 




a=-0.98Jr0.15 

K*-14.8±0.3 

_i I 1 L_ 



-O.a -1 -1.2 

a 



16 



18 



20 



K 

Figure 4. X-band luminosity function of ABCG 209. The solid curve shows the best-fitting Schechter function. The smaller 
panel in the lower-right corner of the figure shows the joint probability contours of K* and a best-fitting values. Black contours 
correspond to la, 2a, and 3a standard probability levels, respectively. The best-fitting value of K* and a are also reported in 
the lower-left corner of the smaller panel. 

4.2 Dynamical properties 

In a previous paper we investigated tlie dynamical status of the galaxy cluster ABCG 209, 
analysing a sample of 119 spectra with 8 A resolution and selecting 112 galaxies belonging to 



the cluster (Mercurio et al. 



2003al ) . Hereafter we will refer to this sample as the "old sample" . 



Among the old sample, 22 galaxies were observed again at higher resolution within the two 
masks together with 29 new galaxies. Combining the old and new redshifts we obtained a 
total sample of 148 galaxies. In Table [1] we report spectroscopic redshift measurements of 
the new observed galaxies. 



4-2.1 Member selection 
In order to select cluster members, we anal ysed the velo city distribution by applying the 



one-d ime nsional adaptive ker nel technique flPisanil Il993l . as implemented by iFadda et al. 



1996 



and 



Girardi et al. 



19961 ). This procedure confirms the existence of a single peak at 



0.209, consisting of 134 cluster members. Moreover, this procedure indicates two other 
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Figure 5. CFHT R-band image of ABCG 209 (North at top and East to left). Galaxies belonging to the cluster are marked 
with crosses, galaxies at 2 ~ 0.257 with circles, galaxies at 2 ~ 0.392 with squares and the remaining isolated galaxies with 
diamonds. 

possible structures at higher redshifts, one at 2; ~ 0.257 with 6 candidate members, and 
another at 2; ~ 0.392 with 4 candidate members. The remaining 4 galaxies are considered 
field galaxies. Although the structure at 2 ~ 0.257 has a high significance (96.6%), the low 
number of galaxies prevents us from assessing the existence of a background cluster. In Fig. [5] 
we show the 134 cluster members as crosses, the 6 galaxies at z ~ 0.257 as circles and the 4 
galaxies at z ~ 0.392 as squares. Figure [6] shows the redshift distribution of the 148 observed 
galax ies. The mean redshift of the cluster as derived by the biweight estimator (Beers et al. 



19901) is < z >= 0.2090 ± 0.0004. 

We further explored the possibility of contamination by interlopers by applying the 



"shifting gapper" algorithm of Fadda et al. (119961 ). This procedure rejects galaxies that are 



too far in velocity-space from the main body of cluster galaxies within a fixed radial bin, 
shif ting a long the distance from the cluster centre. Following the prescriptions of Fadda et 



al. (119961 ) we used a gap of 1000 km s~^ and a bin of 0.6 Xi^q Mpc, or large enough to 
include at least 15 galaxies. In this case eight galaxies were rejected (cf. crosses in Fig. [Tj). 
In order to determine the cluster centre, we applied the two-dimensional adaptive kernel 
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30 



Figure 6. Distribution of redshifts for the 148 observed galaxies. The big and small arrows indicate the mean cluster redshift 
and the redshift of the cD galaxy, respectively. Black areas mark galaxies not belonging to the cluster. 

technique to the galaxy positions. The centre of the most significant peak ( a= 01 31 52.59, 
6= -13 36 45.9) is coincident with the position of the cD galaxy and only 4" distant from 
that of the old sample. 



4-2.2 Velocity dispersion 

We estimated the line-of-sight veloc ity dispersion (LOSVD), a^, using the biweight estimator 
(rostat package; Beers et al. 119901 ). After apply ing the relativistic correction and the usual 
correction for velocity errors (Danese et al. Il980l ). considering only the 126 galaxies selected 
by the "shifting gapper" method, we obtained = 1268^*^84 km s~^, where errors were esti- 
mated with the bootstrap method. If we consider all the 134 galaxies we obtain = I390I9J 
km s~^. This value is slightly larger, but fully consistent with the value computed excluding 
the galaxies identified as outliers by the "shifting gapper" method. 

Velocity anisotropies in the galaxy distribution can create problems in determining the 
LOSVD. For this reason we plot in Fig. [8] the cluster velocity (upper panel) and the velocity 
dispersion (lower panel) vs. projected clustercentric distance. This shows that the a^, com- 
puted by including larger and larger regions, decreases in the central cluster regions, but 
fiattens out in the external regions, suggesting that the value of the LOSVD is no longer 
affected by velocity anisotropies. The rejection of the eight galaxies due to the application 
of the shifting gapper (cf. crosses in Fig. [7]) leads only to a small variation in the estimate 
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Figure 7. Velocity in the cluster rest frame vs. (projected) clustercentric distance for the 134 selected members. The application 
of the "shifting gapper" method with a bin of 0.6 h^Q Mpc rejects the eight galaxies indicated by crosses. 



of cTv and in the velocity disp ersion p rofile (see dot-dashed line in Fig. [8]) so we confirm our 
first result of Mercurio et al. (l2003al ) that the high value of velocity dispersion is connected 
to the strong dynamical evolution of the cluster. In fact, this large value is already reached 
at a radius of 0.4-0.5 hfg^Mpc, where the contamination of interlopers is expected to be 
neghgible. 



4-.2.3 Velocity dispersion of red sequence galaxies 

In order to further investigate the possible enhancement of the velocity dispersion due to 
projection effects from infalling galaxies in the centre of the cluster we obtain the velocity 
dispersion profile considering only galaxies lying on the Colour-Magnitude relation (CMR). 

We obtained the CMR by fitting the photometric data of the spectroscopically confirmed 
cluster members and those galaxies with photometric redshift in the range 0.14 — 0.26 (see 
Sect. 3.2) with a biweight algorithm (Beers et al. 1990), obtaining: 

{R - K)cM = 4.673 - 0.096 -R . (1) 

By using Eq. [H we defined as sequence galaxies the sources lying in the region inside the 
curves: 

(R - i^)seq = iR- K)cM ± 2 ■ ( V^^|+^ + O.OS) , 
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Figure 8. Integrated mean redshift and LOSVD profiles for the 134 galaxies (upper and lower panel, respectively) at a given 
projected distance from the cluster centre, estimated by considering all the galaxies within this radius. The error bars at the 
68% c.l. are shown. In the lower panel, solid and dot-dashed lines give the profile after and before the rejection of eight possible 
interlopers according to the "shifting gapper" method, respectively. The horizontal lines represent X— ray temperature with the 
respective errors transformed in (t„ imposing Papec = 1, where f^spec = (^"i / {kT / ^mp) , with ^ the mean molecular weight and 
nip the proton mass. 



where we take into account the uncertainty at 2a both on the R {(Tr) and o n the K mag nitude 



2002 



20031). 



( ax) as well as the intrinsic dispersion of the CMR (e.g., Merluzzi et al. 

We select 96 spectroscopically confirmed member galaxies lying on the CMR (filled circles 
in Fig. [9]), obtaining a LOSVD for the red sequence galaxies of c'"v,cmr = 

1232^^^ km s-\ By 

selecting galaxies on the CMR, we expect to include only early-type galaxies and to minimise 
the effect of infalling late type galaxies on the line-of-sight velocity dispersion. However the 
measured value of the LOSVD and the velocity profile shown in Fig. [10] are fully consistent 
with those obtained considering all the sample (see dashed line in Fig. [TO!) . 




K 

Figure 9. The R-K/R colour-magnitude diagram of the cluster galaxy population. Open circles indicate galaxies with photo- 
metric redshift in the range 0.15-0.30 and crosses represent the other galaxies with photometric data. The solid line indicates 
the best-fitting C-M relation of Eq. [T] and spectroscopically confirmed member galaxies identified as belonging to the red 
sequence through Eg. 14.2.31 are indicated by solid circles. 




Figure 10. Integrated LOS velocity dispersion profile for the 96 galaxies at a given projected distance from the cluster centre, 
estimated by considering all the galaxies within this radius. The error bars at the 68% c.l. are shown. Solid and dashed lines 
give the profile considering only red sequence galaxies and all the cluster member galaxies after the rejection of eight possible 
interlopers according to the "shifting gapper" method, respectively. The horizontal lines represent X-ray temperature with the 
respective errors transformed in (t„ imposing Pspec = 1, where Pspec = / {kT / ^nip) , with ^ the mean molecular weight and 
rrip the proton mass. 
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4-2.4 Virial radius and mass 

Assuming that galaxies are in dynamical equilibrium with the cluster potential we can deter- 
mine the cluster mass from the knowledge of positions and velocities of the member galaxies 
by applying the virial theorem. Since the application of the virial theorem is meaningful 
only when the system is in dynamical equilibrium within the considered region, a natural 
choice is to compute the cluster mass, Mvij., inside the radius of virialization, -Rvir? within 
which the cluster can be considered to be not far from dynarnical equilibrium. 

Following the prescriptions of Eq. 1 of Girardi & Mezzetti (120011): = 0.17 x a^/H^z), 
where H{z) = Ho(l + z)^/^ (see also Eq. 8 of Carlberg et al. 



19971 for -R200) and is a 



"robust" estimate of the line-of-sight velocity dispersion. We indicate as robust a value that 



is not affected by the presence of velocity anisotropics in t he ga. 



influence the central cluster regions. Following Fadda et al. ll996l ). we check the robustness of 



axy orbits that can strongly 



from the analysis of velocity dispersion profile. Since the profile flattens out in the external 
regions (see Fig. [H]), our estimate of the velocity dispersion is robust (see Sect. 14.2.21) . 

In a ACDM model, with = 0.3, ^]A=0.7 and Hq = 70 km s-^ the value of = 1268^11 
km s~^ lead s to a v alue of the virial radius -Rvir ~ 3.42 h^Q Mpc. In our previous work (Mer- 



curio et al. 



2003al ). we considered a fiat universe with Qm=^- In this case the collapsed 



region for a spherical model corresponds to the region inside Rns- For a comparison with 
our previous estimates we derive i?i78 ~ 2.31 hjQ Mpc. 

On the hypothesis that the galaxy distri butio n traces the mass distribution we obtain 
from the virial theorem (Limber & Mathews Il960l ): 



where v is the galaxy velocity referred to the cluster mean velocity, Vij is the distance between 
any pair of galaxies, a is the global velocity dispersion. 

Assuming that the cluster is spherical, non-rotating system, we can use the projected 
velocity dispersion and the projected virial radius, Rpv- Therefore, Eq. [2] becomes: 

M„. = 1^ . (3) 

In this formula Rpv is the projected radius, equal to twice the (projected) harmonic 
radius: 

Rpv - — ^rr > [V 
where R^^ is the projected distance between two galaxies and is the number of observed 
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galaxies. The value of Rpv depends on size of the sampled region, and on the quality of the 
sampling. In particular, it increases with the cluster aperture A within which it is estimated. 

In our data the cluster is sampled out to a radius of i? ~ l.S/ifo^ Mpc or ~ O.6-R178, but 
we uniformly sample the region Rmax ~ l.O/ifQ^ Mpc or ~ 0.4i?i78. If we consider only the 
region of the cluster that is well sampled, we obtain Rpv = (0.73 ± 0.05)/i7o^ Mpc (the error 
is obtained through the jackknife method) within A = Rmax- Although the cluster is not 
sampled out to Rns, we can estimate Rpv a.t A = Rus theoretically from the knowledge of 
the pa rameters a and Rc of the King-like distribution. We can use Eq. 13 of Girardi et al. 
JlQOsI ). In this case Rpv = Mpc. 

When the system is not entirely included in the observational sample, the usual form of 
the virial theorem 2T + V = should be rep laced by 2T + V = 3PV, where 3PV is the 
surface pressure term (e.g.. The & White Il986l ). Therefore, a surface pressure term C should 
be applied to the mass estimate: 

M, = M„r-C . (5) 

The value of the surface pressure term correction C can be obtained by analysing the 
velocity dispersion profile, but this procedure requires se veral hundreds of member galaxies. 
Combining data of many clusters, Girardi & Mezzetti (120011 ) obtained that velociti es are 
isotropic and that the correction at R ^ R^r is C = 0.2My (e.g. Carlberg et al. 



Girardi et al. 



1997 . 



19981 ). Then, following this approach we can estimate the mass of the system 



inside the collapsed region, assuming a value of 20% for the surface pressure term correction. 
This leads to a value of total mass Mc(< Rns) = ^.U^^qH x 10^^ /i^ Mq, while in a ACDM 
model, with the adopted cosmology Mc(</?vir) = 2.95lo:f8 x 10^^ /i^ Mq- 



5 DYNAMICAL ANALYSIS 
5.1 Velocity distribution 

In order to perform a robust description of the velocity distributi on an d to detect possible 



subclumps we applied the ID KMM algorithm (Ashman et al. Il994j ) to the 126 galaxy 
sample, selected by the "shifting gapper" method, and to the original sample of the 134 
galaxies. The KMM algorithm fits a user-specified number of Gaussian distributions to 
a dataset and assesses the improvement of that fit over a single Gaussian. In addition, 
it provides the maximum-likelihood estimate of the unknown n-mode Gaussians and an 
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assignment of objects into groups. The KMM algorithm is most appropriate in situations 
where theoretical and/or empirical arguments indicate that a Gaussian model is valid, as 
in the case of the cluster velocity distributions, where gravitational interactions drive the 
system toward a relaxed configuration with a Gaussian velocity distribution, but one of the 
major uncertainties of this method is the optimal choice of the number of groups for the 
partition. 

For this reason, we first inves tigate the presence of gaps in the velocity distribution using 



the ROSTAT package (Beers et al. 



19901), which has the advantage of not requiring any a priori 



assumption for the shape of the distribution. Within this package, a gap is defined as the 
difference between two contiguous velocities, weighted by the distances of these velocities 
with respect to their average. We considered only normalised gaps larger than 2.54 since, 
by considering a random sampling of a Gaussian distribution, such gaps appear only in less 
than 1.4% of the cases, independently of the sample size. We found four significant gaps 
indicated by continuous (3.88, 0.05% of the cases), dotted (3.25, 0.10% of the cases), dashed 
(2.59, 1.4% of the cases) and dot-dashed (2.57, 1.4% of the cases) arrows in the lower panel of 
Fig. [m We underline that those gaps are found by using all the 134 galaxies in the sample, 
but identical results are obtained considering only 126 galaxies. 

Using the results of the weighted-gap analysis we divide galaxies into into five groups 
and we apply the ID KMM. From the maximum likelihood statistics, however, we found 
that a mixture of five Gaussians do not give a good description of the data (45.7% c.l.). 
Then we divide galaxies in four groups and we apply the ID KMM, but from the maximum 
likelihood statistics, we found that the data can be described as a mixture of four Gaussians 
at only 89.7% c.l. On the other hand, from the maximum likelihood statistics we found 
that a mixture of three Gaussians is the best description of the velocity distribution (at 
98.4% C.L). In the upper panel of Fig. [11] we plot the velocity distribution of galaxies with, 
superimposed, the three Gaussians corresponding to the three identified clumps. The most 
prominent gap detected by the rostat package is between the clumps 1 and 2 identified by 
the KMM algorithm. Using the one dimensional KMM algorithm we assigned the member 
galaxies to individual groups (ni=16, n2=102, and n3=16 members at mean redshift zi = 



-47 

-58 



0.1988, Z2 = 0.2090, and Z3 = 0.2172) and we estimate the velocity dispersion of o"vi = 323" 
km s~^, a^2 = 847^49 km s~^, and a^s = 289llll km s~^(see FigJTTl) for each clump. We note 
that the old sample of 112 member galaxies did not allow us to associate galaxies to the less 
populated substructures and then to obtain a reliable estimation of their velocity dispersions. 



Global Properties of ABCG 209 25 



15 - 




5.8x10* 6x10* O.axlO* 6.4x10* 6.6x10* 

V [km s 'J 



Figure 11. Distribution of radial velocities for member galaxies. Upper panel: velocity histogram with a bin of 250 km s~ , 
with the three solid Gaussians showing the best separation into three clumps as identified with the one-dimensional KMM 
test (see text), and the dotted Gaussian indicating the best-fit single Gaussian function to the velocity distribution. Lower 
panel: stripe density plot . The arrows indicate the positions of the three significant gaps in the velocity distribution of member 
galaxies (see Beers et al. ll99ll V 



Figure fT2l shows the spatial distributions of the three clumps. Differently from the results 
of our first paper, we find no evidence for spatial segregation for the three clumps, accord- 
ing to the two-dimensional Kolmogorov-Smirnov test (hereafter 2DKS-test; cf. Fasano & 



Franceschini 



19871 . as implemented by Press et al. 



1992). 



5.2 3D substructures 

In order to check for the presence of substructure, we combined velocity and p osition in 



formation by computing the A-statistic devised by Dressier & Schectman (119881 ). This test 
is sensitive to spatially compact subsystems having an average velocity and/or a velocity 
dispersion different from the global mean quantities. We found a value of 198 for the A 
parameter, which gives the cumulative deviation of the local kinematical parameters (ve- 
locity mean and velocity dispersion) from the global cluster parameters. The significance 
of substructure was checked by running 1000 Monte Carlo simulations, randomly shuffling 
the galaxy velocities, obtaining a significance level of 99.8%. This indicates that the cluster 
has a complex structure. In Fig. [13] we plot the member galaxies marked by circles whose 
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Figure 12. Distribution of member galaxies separated into three clumps aceording the one-dimensional KMM test. The plot 
is centred on the cluster centre. Open circles, dots and crosses indicate clumps 1, 2 and 3, having low, intermediate and high 
velocities, respectively. 



diameter is proportional to the deviation 5 of the individual parameters (position and ve- 
locity) from the mean cluster parameters. We confirm with this larger sample the results of 
the previous paper. In particular, a group of galaxies with high velocity in the external East 
cluster region and another group near the cluster centre could cause the large values of S. 
We also underline that we obtain a value of 199 for the A parameter with a significance level 
of 99.9%, by considering only the 126 galaxies selected with the shifting gapper method. 

Since the technique by Dressier & Schectman does not allow a direct identification of 
galaxies belonging to the detected substructure, we apply the three dimensional version 
of the KMM test using simultaneously galaxy positions and velocities. Wc use the galaxy 
separation obtained in the ID analysis as first guess of the 3D analysis. The algorithm fits a 
three-group partition at 98.6% c.L, with a galaxy distribution similar to those obtained by 
the ID analysis (ni=16, n2=108, and n3=10 members). On the other hand the 3D KMM fits 
a four-group partition at 100.0% c.L, with a partition similar to those obtained in the gap 
analysis performed with the rostat package (ni=16, n2=32, n3=29, and n4=47 members). 
However, this test seems to be too sensitive to changes in the initial conditions, in fact with 
a particular choice of the initial group partition it is possible to obtain a separation into two 
groups significant at 98.0% c.L 
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Figure 13. Spatial distribution on the sky of the 134 cluster members, each marked by a circle: the larger the circle, the 
larger is the deviation 5 of the local parameters from the global cluster parameters. The figure shows evidence for substructure 
according to the Dressier & Schectman test. The plot is centred on the cluster centre. 



6 STRUCTURAL PROPERTIES OF EARLY-TYPE GALAXIES 
6.1 Structural parameters 

We derived the structural parameters, namely half-light radius, re, average surface brightness 
inside r^, /ig, and Sersic shape parameter, n, in the R and K band for a large sample of 
galaxies belonging to ABCG 209. 

Structural parameters were deriv ed fo r galaxies in the R-, and K-band, by using the 



2DPH0T package (La Barbera et al. l2008|) 



We considered the sample of spectroscopically confirmed cluster members, complemented 
with the galaxies having photometric redshift in the range of 0.14 — 0.26 (see Sect. 13. 2p . We 
found that for Structural parameters were derived for galaxies with R < 21.0 mag and 
K < 17.5, i.e. the limits of spectroscopic observations, leading to a total sample of 327 
galaxies. Numerical simulation were performed in order to check the reliability of the derived 
structural parameters, finding that the typical errors were (5(logre) ~ 0.14 and 6{< /io >) ~ 
0.55 on the effective radius and surface brightness respectively. 

Structural parameters were derived by fitting galaxy images with seeing-convolved Sersic 
models. The point-spread functions (PSFs) were derived locally from the stars in the R and 
K images and were modelled by a sum of 2D Moffat functions, also taking into account 
deviations of stellar isophotes from circular symmetry. 
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Details on the dif ferent a spects of the der ivatio n of s tructural parameters may be found 
in La Barbera et al. (20o3), (2003a), (2003b) and (2008). 

The structural parameters are given in Table [3] for the whole sample of 327 galaxies with 
photometric and spectroscopic redshifts. We remark that this constitutes so far the largest 
sample of (optical plus NIR) structural parameters of galaxies belonging to one cluster. 



6.2 Kromendy relation and Photometric Plane 

ETGs define a remarkable correlation between the effective radiu s r^ a nd the mean surface 
brightness < fie>, known as Kormendy relation (KR, Kormendy 119771 ) : 

< /ie > = a + (3 logre . 

In order to study the Kormendy relation (hereafter KR), we selected the population of 
spheroids on the basis of the shape of the light profile, as parameterised by the Sersic index 



n. We classified as spheroids the galaxies with n > 2 , corresponding to o 



fraction greater than ~ 20% (see Saglia et al. 



jjects within a bulge 



1997: van Dokkum et al. 



1998f ). Moreover, 17 



and 54 objects in R and K-band respectively, with small radii were excluded by selecting 
only galaxies with r^. greater than 1 pixel. 

Again, we point out that the present samples constitute the largest data set of galaxies 
belonging to one cluster for which the KR is obtained. In fact 240 and 227 galaxies are iden- 
tified as spheroids in R and K band respectively, with spectroscopic redshift or photometric 
redshift in the range 0.14 — 0.26 (see above). 

Si nce se lection effects can strongly affect the estimate of the KR coefficients (see Ziegler 



et al. 



19991 ). we fitted the logre-< /ie > sequen ces by introducing a modified least-squares 



(MLS) procedure (see La Barbera et al. l2003bl ). which corrects the bias due to the different 
completeness cuts in magnitude of each sample. The fi tting coefficients were derived by 
applying the bisector regression (see Akritas & Bershady Il996l ). 

Figure [H] shows the KR for R (left panel) and K band (right panel), reporting also the 
correlation between (51ogre and 5< fie >, with the typical values of the uncertainties on the 
effective parameters. We obtain zero points a = 19.19 ± 0.07 and a = 15.68 ± 0.07, and 
slopes P = 3.07 ± 0.16 and (3 = 3.26 ± 0.19 for R and K bands, respectively. 

The difference in the zero points is fully consistent with the colour R-K of galaxies, 
according to the CM relation (Fig. [9]). 

The slope of the KR in R band is fully consistent with that derived in La Barbera et al. 
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Table 3. Structural parameters of galaxies in ABCG 209, according to the photometric redshifts. Right ascension and declina- 
tion (Col. 1, Col. 2) effective radius, mean surface brightness and Sersic index for R (Col. 3, Col. 4, Col. 5) and K band (Col. 6, 
Col. 7, Col. 98) 
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Figure 14. Kormendy relations for ABCG209 in R (right panel) and K band. MLSiogr^ is represented by long dashed lines. 
The short dashed line indicate the cuts in total magnitude. The correlation of the uncertainties on logrc and < /.tc > is shown 
by the ellipses (1 o" confidence contours) in the lower left of each panel. 



fl2003bl ) for a smaller (N=81) and shallower (R<20.1) sample of galaxies in ABGC209, thus 
confirming that the R-band KR is invariant in the redshift range from z=0.023 to z=0.64. 
For what the K-ban d KR is concerned, our slope is fully consistent with that found by La 



Barbera et al. (l2004l ) for the galaxies in the cluster A2163B at z~0.2. 



The structural properties of early-type galaxies define a plan e in the three-dimensional 



space (logre,logn, < /le >), the so-called PHP (see Graham |2002| and references therein): 



logrg = a logn + b < yUg > + c . 

We derived coefficients of the P HP by using the corrected orthogonal weighted least- 
squares fit (La Barbera et al. |2005| ). This method, treating equally all the variables, gives 
much more robust estimates for the coefficients of the PHP and is less sensitive than others 
to selection effects. 

The distribution of galaxies in the space of structural parameters is shown in Figs. [15] and 
[T6l for R and K band respectively. We show the distribution of galaxies in the logr^.— < yUg >, 
logn— < fie > and logn — logve planes, and an edge projection of the plane. Galaxies follow a 
well-defined PHP at z= 0.2, with Sersic indices that increase towards lower surface brightness 
values and larger effective radii. 

In Table H] we report the coefficient of the PHP for R and K band. The coefficients of 
the R- and K-band PHP are consistent (at lo") with those obtained by La Barbera et al. 
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Table 4. Coefficients of the PHP. 



Band a b c (Tlo^ 



R 
K 


0.964 ±0.130 
0.711 ±0.098 


0.209 ±0.020 
0.227 ±0.022 


-4.376 ± 0.401 
-3.871 ±0.404 


0.162 ± 0.014 
0.167 ±0.016 
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Figure 15. Photometric Plane of ABCG 209 in the R band. The upper panels and the lower-left panel show the logre— < /ie >, 
logn— < fie > and logn — logr^ projections of the PHP, respectively. The lower-right panel shows the edge-on view of the 
R-band Photometric Plane. The plots show the N= 240 galaxies of the R-band sample. 



(120051 ) for the cluster MS 1008 at z~0.3 and by Graham et al. (120051 ) for nearby clusters, 
confirming that PHP seems to be independent of redshift. 



6.3 The intrinsic dispersion of the PHP 

In Table m we showed that the PHP has a sma ll intrinsic dispersion {(Jiogre = 0-16 — 0.17) 



both in R and K band. La Barbera et al. 



20051 pointed out that stellar populations can be 



the origin of the dispersion about the plane. However they have shown that the residuals 
of the I-K versus K colour-magnitude relation do not correlate with the residuals about the 
PHP. 



32 A. Mercurio et al. 




-0.5 



0.5 1 
Log (kpc) 



1 .5 



1.5 



Log n 




Log n 



1 .5 



^0.5 



0.5 1 

Log r (Kpc) 



1 .5 



Figure 16. Photometric Plane of ABCG 209 in the K band. The upper panels and the lower— left panel show the logre— < lie >, 
logn— < fie > and logn — logre projections of the PHP, respectively. The lower-right panel shows the edge-on view of the 
K-band Photometric Plane. The plots show the N= 227 galaxies of the K-band sample. 



The dispersion about the PHP could be due to the combined effect of age, metalhcity, 
and also a/Fe enhancement, that cannot be distinguished only with the colour-magnitude 
relation. For this reason it is crucial to use line strength to address this issue. 

In order to distinguish between the contribution of age, metalhcity and a/Fe enhance- 
ment, we compared the line indices obtained in Mercurio et al. ( 120041 ) with the residuals of 
PHP. We analyse a subsample of 83 and 91 ETGs, in R and K band respectively , for w hich 
we have these measurements. In particular we use H/j, [MgFe]' (Thomas et al. Il986l ) and 
Mgb/<Fe>, as the b est in dicators of age, metalhcity and a/Fe enhancement, respectively 



(see also Puzia at al. 



20051). 



Figure [T7| shows that the H^ (upper panel) and [MgFe]' indices (central panel) do not 
correlate with either the R-band PHP (filled circles) or the K-band PHP (open circles). The 
dispersion of the PHP seems instead related to Mgb/<Fe> measurements (lower panel). 

We perform a least-square fit on the data, obtaining -0.051±0.461 and -0.261±0.705 for 
the age and metalhcity correlation slopes, while for the a/Fe we obtain 0.702±0.258. 
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Figure 17. Correlation between PHP residuals and equivalent width of the (upper panel), strength of [MgFe]' (central 
panel) and Mgb/<Fe> (lower panel) for R (filled circles) and K (open circles) band. Continuous lines represent the least-square 
fit on the data. 



This suggest that the chemical evolution of galaxies could be responsible for the intrinsic 
dispersion of the PHP, as argued by La Barbera et al. (120051 ). 



7 SUMMARY AND CONCLUSIONS 

In order to study the internal dynamics of the rich galaxy cluster ABCG 209, we have 
analysed spectra for a total sample of 148 galaxies in the cluster region, of which 134 are 
candidate cluster members having z~0.209. In a previous paper we investigated the dynam- 
ical st atus of this cluster, selecting 112 galaxies belonging to the cluster (Mercurio et al. 
2003al ). We complement the previously described dataset with medium resolution spectra 
for 29 new cluster galaxies. We also re-observe 22 cluster members contained in the already 
available data sets. The enlarged data sets allow us to better distinguish among galaxies 
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belonging to different substructures, to derive their individual velocity distributions and to 
study possible spatial segregation. 

We extend the existing optical (i?, V and R bands) photometry with i^-band observa- 
tions which allow us to estimate reliable photometric redshifts {5z/l + z ^ 0.07) for 399 
galaxies brighter than R=21 mag in the cluster field and to derive the NIR LF. 

The Schechter function provides a good description of cluster galaxy counts, with a = 
—0.98 ± 0.15 and K* = 14.80 ± 0.30 mag, values consistent with previous studies of NIR 



luminosity functions of clust er galaxies at z~0.2 (Ide Propris et al. 



2000 



La Barbera et al.l 



200^ . The only magnitude bin showing a larger deviation is that at 



1998 



Andreon fc Pello 



K = 17.5 mag, where a marginal (~ 1.5cr significant) deficiency of cluster galaxies is found, 
that corresponds to the magnitude bin (20 < i? < 21 mag) where Mercurio et al. (j2003bl ) 
find some indication of a dip in the optical luminosity functions of ABCG209. 
ABCG209 is characterised by a very high value of the LOSVD: 



1268tll km s-\ 



that results in a virial mass of M„ 



2.95 



+0.80 
'0.78 



X lO^^/iyn hyn Mm within R, 







3.42/ifo^ 



Mpc. The analysis of the velocity dispersion profile shows that such a high value of is 
already reached in the central cluster region suggesting the possibility that a mixture of 
clumps at different mean velocities causes the high value of the velocity dispersion and the 
virial mass. In fact, we found that a mixture of three Gaussians is the best description of the 
velocity distribution (at 98.4% c.l.). Assigning the member galaxies to individual groups we 
found the main clump (ngaj=102) at mean redshift z = 0.2090 with a velocity dispersion of 

= 847^49 km s~^and other two groups of 16 galaxies each at z = 0.1988 and z = 0.2172 
with (Tv = 323^58 km s^-'^and = 289^54 km s~^respectively. 

This observational scenario confirms that ABCG 209 is presently undergoing strong dy- 
namical evolution with the merging of two or more subclumps, in agreement with the recent 
detection of a radio halo (Giovannini et al. |2006| ). On the other hand, there a re dis crep- 
ancies betwee n the dynamical. X-ray and lensing analyses. While Smith et al. (120051 ) and 
Zhang et al. (120071 ) found no evidence for multimodality in the dark matter distribution 
and subclumps in the XMM data respectively, an irregular X- ray mo rphology was identified 



by Rizza et al. ( 



Smith et al. 



2005 



998l ) in ROSAT data and by Mercurio et al. l2003al (see their Fig. 13) and 



see their Fig. 6) in Chandra observations. In particular, the Mercurio et al. 
(I2003al ) analysis recognised two significant substructures. The principal one centred on the 
cD galaxy and the second about 50" East of the cD, t his secondary clump being coincident 
with the Eastern clump detected by Rizza et al. ( 19981 ). Since these two detected clumps are 



Global Properties of ABCG 209 35 

close to each other, probably the absence of evident subclumps in the XMM-Newton data 
is due to the limited spatial resolution of this instrument. 



Smith et al. (120051 ) suggest that this discrepancy could be related to time evolution and 
that s ufficie nt time has elapsed after the merger. However, as also pointed out by Smith 
et al. (120051 ). numerical simulations suggest that both gas dynamics and substructure in 
the dark matter distribution may persist as long as ~5 Gyr after a cluster merger, and 
that the relaxation time for these two matter c ompo nents may be comparable. Moreover, 
the detection of a radio halo (Giovannini et al. |2006| ) supports a recent merging event. In 
fact numerical simulation indicate that the mergers could generate strong fluid turbulence, 
supplying energy to the electrons, which then radiate in the interstellar r nedium of radio 
halos, but the time during which the process is effective is ~10^ yr (Feretti 120071 ). 

In independent lensin g ana. 



yses of the optical CFHT images, Paulin-Henriksson et al. 
(j2007l ) and Bardeau et al. (j2007l ) measured virial masses of M200 = ^.Tt^^r x 10^^ hfg^ and 
M200 = 7.2 ± 2.0 X 10^^ hfo^ M0 and scaled velocity dispersions = 924 ± 84 km s"^ and 
av = 813 ± 70 km s~^ respectively, which are consistent with the value obtained for the 
main clump from the dynamical analysis. For this reason we suggest that the observed weak 
lensing masses are associated with the dark matter halo of the most prominent clump of 
ABCG 209. The clumps are not spatially segregated so they would not be separable in the 
dark matter maps which measure only the projected mass distribution, and the combined 
contribution of the secondary clumps is expected to be much lower than that from the 
primary one. 

We have derived structural parameters of galaxies in the R and K bands in order to 
obtain the Kormendy relation and the Photometric Plane for ABCG 209, based on a total 
sample of 240 and 227 spheroids in optical and NIR, respectively. Spheroids deflne a tight 
sequence in the plane of the effective parameters with a slope /5 ~ 3 in bot h bands . This 
slope of the KR is fully consistent with that derived in La Barbera et al. (j2003bl ) for a 
smaller (N=81) and shallower (R<20.1) sample of galaxies in ABGC209, thus conflrming 
that the R-band KR is invariant in the redshift range from z=0.023 to z=0.64. For what the 
■ "C-ban d KR is concerned, our slope is fully consistent with that found by La Barbera et al. 



20041 ) for the galaxies in the cluster A2163B at z~0.2. 



The cluster ETGs at follow a tight correlation among logrg, < fi >e and log n, 

with an intrinsic dispersion of ~0.17 dex in both optical and NIR wavebands. This dispersion 
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is fully consistent w ith th at found by La Barbera et al. (120051 ) for the cluster MS1008 at 
2:^0.3 and Graham (120021 ) for ETGs in nearby clusters. 

In order to investigate the origin of the intrinsic scatter of the PHP we have analysed 
the PHP residu als ve rsus line-strength indices. We compare the line indices obtained in 
Mercurio et al. (120041 ) with the PHP residuals of member galaxies for a subsample of 83 
an d 91 E TGs, in R and K band respectively. In particular we use H^, [MgFe]' (Thomas et 



al. 



19861 ) and Mgb/<Fe>, a s the best indicator s of a ge, metallicity and a/Ye enhancement, 



respectively (Thomas et al. 



1986 



20051 ) . The PHP residuals do not correlate 



Puzia at al. 

with age and metallicity, while there is a correlation with a /Ye enhancement (Fig. [T71) . This 
could imply that the scatter of the PHP is due to variations in stellar population parameters, 
in particular it could be due to variations in the chemical evolution of early-type galaxies 
brought about by their merging histories. 
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